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Dear Sirs, 

I am contacting you with respect to a patent application. No: 10/813,499 entitled 
"Automatically conforming intraoral dental radiographic sensor". Fundamental parts of the 
features claimed in this applications are not novel, as you may see from two publications 
included with this inail. More specifically, besides a Gernian Patent application in 2001 
(which, of course is not relevant for the US) the correction of projective distortion using 
fiducial spherical reference bodies has been described in detail in an internationally published 
paper in 2002 [1]. We even extended the method for complete registration of projection 
geometry in another paper published in a high-standard (US-based) Journal in 2004 [2]. 
I believe, these papers should be considered in the patenting process for this particular 
application, particularly when novelty is targeted. 

For your information, I enclose both papers together with the patent application with this 
mail. 

Thank you very much for your kind attention to this objection. I would be very grateful, if you 
could send me some information on the final decision on that issue. 



Sincerely, 




1. Schuize, R., d'Hoedt, B., (2002). A method to calculate angular disparities between object and receptor 
in "paralleling technique". Dentomaxillofac Radiol; 31: 32-38. 

2. Schuize, R., BrUllmann, D. D., Roder, F., d'Hoedt, B., (2004). Determination of projection geometry 
from quantitative assessment of the distortion of spherical references in single - view projection 
radiography. Med Phys; 3 1 : 2849-2854. 
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ABSTRACT 



A system is described for obtaining intraoral radiographic 
images using an anatomically conforming sensor and a 
method is described for correcting the projective distortions 
caused by the inclination of the plane of the imaging sensor 
relative to the x-ray beam. The purpose of the described 
sensor system is to reduce patient discomfort as compared to 
conventional planar intra-oral x-ray sensors, which tend to 
impinge on sensitive tissues within the oral cavity. An 
additional purpose of the described system is lo allow for 
better coverage of the subject by allowing the sensor to fit 
closer against the surrounding tissues. 
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ANATOMICALLY CONFORMING INTRAORAL 
DENTAL RADIOGRAPHIC SENSOR 

FIELD AND BACKGROUND OF THE 
INVENTION 

[0001] The present invention relates4o=^n intraoral sensor 
design to be applicable to radiographic imaging sensors such 
as fihn, storage phosplror-dwices or solid-stale direct cap- 
lure devices (e,g. CGE?s or CMOS based devices) or any 
other imaging sensor that can be used to produce radio- 
graphic images. Such sensors may be components of radio- 
graphic imaging systems that may include computers and 
video monitors for displaying radiographs. Tlie present 
invention is also directed toward a method of correcting 
distortions appearing/in the images^of such radiographic / 
imaging systems caused by the non-perpendicularily of the ' 
x-ray source relative to the imaging surface. Such analyses 
may be performed on a PC, which in turn may be part of a 
digital radiographic imaging system. The illuminating x-ray 
radiation source for such sensors and radiographic imaging 
systems is not the subject of the present invention. 

[0002] Radiographs are fundamental to most dental diag- 
nostic procedures. However, a common complaint and prob- 
lem during radiographic exams is patient discomfort during 
the placement of radiographic sensors within the mouth. The 
majority of these complaints involve the placement of the 
radiographic sensor in the posterior maxillary and mandibu- 
lar arches of the patient. This problem is primarily due to the 
limited space available for proper placement of the sensors 
within these regions. This has been a problem since the 
inception of dental radiography using standard x-ray fihn 
technology. 

[0003] A common method used to reduce such patient 
discomfort is to bend portions of the film package to 
conform to the anatomic shape of the mandibular or max- 
illary arches. However this bending contributes to the prob- 
lems of image distortion of the intended subject as described 
below. 

[0004] Recently, solid-state x-ray sensors have been 
developed that replace film. The patient discomfort problem 
for these sensors is even greater because these devices are 
rigid by nature and cannot be bent like film. With these rigid 
sensors, the only recourse is to place the sensor at a more 
comfortable angle and position for the patient, but such 
positioning is normally not optimal for a diagnostic radio- 
graph of the intended subject. One study has shown that 
using such rigid devices results in more sensor placement 
errors than conventional film, thus requiring more.retaHes to 
producecUnically acceptable radiographs (C.H. Versteeg, el. 
al., Art evaluation of periapical radiology with a charge- 
coupled device; Dentomaxillofacial Radiology, 1998, 27: 
97-101). 

[0005] Both of these approaches, using film and rigid 
sensors, result in a distorted radiographic image of the 
intended subject when the x-ray beam is not perpendicular 
to the imaging plane. 

[0006] Thus, a need exists to provide for a more comfort- 
able sensor configuration while correcting for the projective 
distortions of the images produced when the x-ray beam is 
not^perpendicular to the imaging plane. " 



OBJECTS AND SUMMARY OF THE PRESENT 
INVENTION 

[0007] Accordingly, an object of the present invention is to 
provide a configuration for a radiographic sensor that allows 
for increased patient comfort. 

[0008] Another object of the present invention is to pro- 
vide a configuration for a radiographic sensor that allows for 
improved coverage of the radiographic subject, 

[0009] Yet another object of the present invention is to 
- provide for the incorporation of fiduciary elements within a 
radiographic sensor that allows for the measurement of the , 
magnitude of image distortion produced when the imaging- 
x-ray beam is not perpendicular to the sensor surface. ' 

[0010] Finally, another object of the present invention is to 
provide a method for the correction of distortions in a 
radiograjjliie image that results when the x-ray beam is not 
perpendiciiiar'tb the plane of the sensor. 

[0011] The present invention is embodied in a radiation 
sensor that has a specific geometric shape or configuration 
that not only provides more comfort to a patient while in 
intraoral use, but one that is capable of providing a more 
accurate image of the area of concern by eliminating dis- 
tortions of the resultant image normally found in other 
modem radiographic imaging systems and sensors. Such 
systems and sensors typically operate in conjunction with 
personal computer systems (PC's) whose function is to serve 
as a platform to receive, analyze, display and archive image 
data received from the sensor systems. 

[0012] The preferred geometric configuration of a x-ray 
radiation imaging sensor that embodies the present invention 
is unique among known intraoral radiographic imaging 
sensors in that it divides the radiographic imaging sensor 
into several discrete sections, each section imaging a dif- 
ferent region of the subject. 

[0013] Furthermore such sections are arranged in a fashion 
that conforms to the anatomic curvatures of the human 
maxilla and mandible while maintaining a diagnostically 
useful image and contiguous coverage of the subject. To 
accomplish this the different sensor sections may be oriented 
in different directions, resulting in a non-planar shape. 

[0014] Also, the sections of the sensor are abutted together 
in a fashion that yields contiguous coverage of the subject 
area. Finally, each section has incorporated within it radio- 
opaque fiduciary elements of known shape and size whose 
projected radiographic images indicate the orientation of the 
x-ray beam relative to the surface of each of the sensor 
sections. 

[0015] Another embodiment of the present invention 
involves a method of processing of image data produced by 
the hardware that also embodies the present invention. This 
processing method includes-steps that use the image of a 

rfidiiciary element, whose shape and size is known, to correct . 

Uhe projective distortions in the resultant radiographic image 
caused by non-perpendicularity of the x-ray beam relative to 
the surface of the sensor section. 

[0016] Such inclination of the incident x-ray beam causes 
a projective distortion of the radiographic image of the 
fiduciary element, as well as the radiographic image of the 
intended subject. Because the shape and size of the fiduciary 
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element is known, the precise orientation of the x-ray beam 
._can be derived from the nature of the distorted image using 
principles of projective geometry. In turn, the resultant 
distorted radiographic image can be restored to ils undis- 
torted state from Imowledge of the beam orientation through 
_a reverse projective transformation of the image, again using 
principles of projective geometry. This method of the inven- 
~ tion is concerned with image distortion and restoration and 
is independent of the number of sensor sections and applies 
even to sensors consisting of only one section. 

[0017] In general, one embodiment of the present inven- 
tion is found in a radiation sensor that is intended for use 
with a digital intraoral radiographic imaging system for 
placement in the mouth of a patient for production of 
radiographs of teeth and supporting structures. Such a sensor 
includes a housing containing at least two generally planar, 
radiation detectors abutting at a non-zero angle forming a 
generally contiguous convex or concave image capturing 
siu-face oriented toward a radiation source, The housing 
shape is selected to conform to the anatomical curvatures of 
the human maxillary and mandibular arches of the average 
patient. At least one radio-opaque fiduciary elecnent of ' 
known shape, size and location is embedded on, in or just 
under the housing intermediate to the radiation source and to 
the surface of each of the radiation detectors. The fiduciary 
elements cast a projected image on the surfaces of the - 
radiation detectors when illuminated by an x-ray source^ for 
u^ in eliminating distortions in the image captured by the 
sensor. 

[0018] Another embodiment of a sensor containing the 
present invention includes a radiation sensor also intended 
for use with a digital intraoral radiographic imaging system 
for placement in the mouth of a patient for production of 
radiographs of teeth and supporting structures that conforms 
to the anatomic oral structure of the patient in a manner that 
provides less distortion in the detected image than that of 
current sensors. In this embodiment the sensor is enclosed in 
a housing containing at least two generally planar, radiation 
detectors abutting at a non-zero angle to form a faceted, 
generally contiguous convex or concave image capturing 
surface oriented toward a radiation source. I'he angle is 
selected such that the sensor housing conforms to the 
anatomical curvatures of the human maxillary and mandibu- 
lar arches of the average patient. In this embodiment, each 
adjoining pair of generally planar radiation detectors are 
held in a fixed angular relationship to one another or flexibly 
joined so that the angle at which they abut can be changed 
to allow the housing to conform more easily to the anatomi- 
cal curvatures of the human maxillary and mandibular 
arches of the average patient. 

[0019] Additionally, the present invention is also found in 
a method for correcting distortions in a radiation sensor used 
with a digital intraoral radiographic imaging system for 
placement in the mouth of a patient for the production of 
radiographs of teeth and supporting tissues. A method 
embodying the present invention includes the steps of: 
[0020] Providing a housing containing at least two gener- 
ally planar, radiation detectors abutting at a non-zero angle 
to form a faceted, generally contiguous convex or concave 
imaging surface oriented toward a radiation source, the 
angle being selected such that the sensor housing conforms 
more easily to the anatomical curvatures of the human 
maxillary and mandibular arches of the average patient; 
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[0021] Placing at least one radio-opaque fiduciary element 
of known shape, size and location intermediate to the 
radiation source and to the surface of at least one of the 
^radiation detectors; 

[0022] Exposing the surfaces of the radiation detectors and 
the fiduciary element to the radiation source to project an 
image of the fiduciary element and of the teeth and support- 
ing structures onto the radiation detectors; 

[0023] Capturing and digitizing for analysis the sensor 
output representing the projected image of the fiduciary 
elements and the teeth and supporting structures detected by 
the radiation detectors; 

[0024]. Analyzing the digitized data representing the 
image of the fiduciary element to determine the magnitude 
of distortion of the projected fiduciary element's shape onto 
the surfaces of the > radiation detectors, due to the non- 
perpendicularity of the radiation source relative to the sur- 
faces of the radiation detectors in comparison to that of an , 
ideal undistorted image of the fiduciary element projected 
, onto the surfaces of the radiation detectors by expo.sure to a 
radiation source that is perpendicular to the surfaces of the 
radiation detectors, and determining a corrective transfor- 
mation that transforms the distorted image of the projected 
fiduciary element into that of an undistorted image^of the 
fiduciary elements; and, ~ 

[0025] Applying the corrective transformation to the digi- 
tized image data to transform the distorted projected image 
of the teeth and supporting stmctures to that of an undis- 
torted projected image of the teeth and supporting strucmres. 

[0026] Various additional objects and advantages of the 
present invention will become apparent from the following 
detailed description. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0027] FIG, Ifl is a cross-se<Aional view through the 
posterior portion of the maxilla illustrating an anatomically 
shaped radiographic sensor embodying the present inven- 
tion; 

[0028] FIG. li? is a cross-sectional view through the 
posterior portion of the maxilla showing the placement of a 
conventional radiographic sensor and illustrating the poten- 
tially painful impingement of the sensor on the soft tissues 
of the palate; 

[0029] FIG. 2fl is a cross-sectional view through the 
posterior portion of the maxilla and an internal cross- 
sectional view of a section of an intraoral sensor illustrating 
an alternate anatomically shaped radiographic sensor with 
housing, radiation detectors and incorporated fiduciary ele- 
ments embodying the present invention; 
[0030] FIG. 26 is a frontal view of an alternate anatomi- 
cally shaped radiographic sensor similar to that shown in 
FIG. 2j illustrating a possible spatial distribution of fidu- 
ciary elements within a sensor embodying the present inven- 
tion; 

[0031] FIG. 3fl is an idealized schematic showing the 
projective geometrical relationship of a fiduciary clement 
onto a plane; 

[0032] FIG. 36 is a schematic showing the internal geo- 
racGic arrangeinent of a fiduciary element with respect to the 
suffaoe'oftte radiation detector in a section of a sensor and 
its projected image thereon; 
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[0033] FIG. 4fl is an image of a spherical fiduciary ele- 
ment projected onto a sensor's radiation detector's surface; 

[0034] FIG. 4b is an image of the spherical fiduciary 
element projected onto a sensor's radiation detector's sur- 
face as in FIG. 4a after the application of an edge filtering 
transformation to the image data; 

[0035] FIG. 4c is the image of FIG. 4a after the image has 
been subjected to threshold and thinning filters; 

[0036] FIG. Sa illustrates a distorted image of maxillary 
teeth due to excessive x-ray beam inclination; and, 

[0037] FIG. Sb illustrates a corrected image of the dis- 
torted image of maxillary teeth as shown in FIG. 5fl using 
the distortion of the projected images of the spherical 
fiduciary elements as a guide in correcting the distortion in 
the projected image plane. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENT OF THE INVENTION 

[0038] Referring now to the Figures, a sensor embodying 
the present invention is shown generally at 100 in FIG. la. 
The simplest implementation of an anatomically shaped 
radiographic sensor 102 embodying the present invention 
has two abutted imaging planes 104, 106, respectively, as 
shown in FIG. la which illustrates a cross-sectional view 
through the posterior portion of a patient's maxillary arch. A 
sensor holder/positioner 108 is shown for retention between 
a patient's teeth to hold sensor 102 in position during 
imaging procedures. 

[0039] The specific angle of the two imaging planes 104, 
106 relative to each other is not critical but empirical testing , 
indicates that angles of about 20-40 degrees result in sig- 
nificant patient comfort relative to a single flat plane of 
equivalent imaging area. This configuration has been shown 
to significantly improve comfort and improve subject cov- 
erage compared to a standard sensor. 

[0040] FIG. 1/? is a similar view of a conventional radio- 
graphic sensor 110 and its corresponding sensor holder/ 
positioner 112 illustrating the potentially painful impinge- 
ment of the sensor on the soft tissue of the maxillary palate 
and the unfavorable imaging placement due to the increased 
distance of the sensor relative to the teeth and supporting 
structures. 

[0041] FIG. 2ii shows the placement of fiduciary elements 
114 relative to the surfaces of the radiation detectors 116 of 
the sensor 102 and its housing 120. Criteria for such 
elements include: 1) a size sufficiently small to minimi/x 
obstruction of critical features in the radiograph, 2) a size 
large enough to allow precise length measurement of the 
element relative to the resolution of the sensor, 3) radio- 
opacity of the element and, 4) as a manufacturing consid- 
eration the fiduciary element should be symmetrical in shape 
so that its precise orientation on the sensor is not critical. 

[0042] FIG. 2b shows a frontal view of the preferred 
embodiment, illustrating the placement of the fiduciary 
elements 114 relative to the sensor surface 118. Size and 
placement of such fiduciary elements 114 should strive to 
ensure imaging of at least one element per section of sensor, 
despite misaiming of the illuminating x-ray beam while 
minimizing obscuration of the intended subject. 



[0043] Thus, in the preferred embodiment a 1 mm diam- 
eter stainless steel symmetric sphere has been selected as the 
fiduciary element 114 being illustrated. This is because 
stainless steel offers sufiScient radio-opacity to the imaging 
X-ray beam, and also because currently available high- 
resolution CCD sensors typically have a pixel resolution of 
0.022 millimeters (22 microns). Hence, the size measure- 
ment of the element should be accurate to about two percent 
(more than adequate for most dental diagnostic applica- 
tions). Film offers even higher measurement accuracy. 

[0044] The image distortions caused by non-perpendicu- 
larity of the illuminating x-ray beam relative to the surface 
of the radiation detector belong to a general class of math- 
ematical_functions known as projective transforms. More- 
over, if the illuminating source is relatively- far from the 
intended subject (far field approximation), a good approxi- 
mation for dental imaging, then the distortion belong to a 
^special sub-class of projective transforms known as afi&iie/ 
^transforms. In the case of the preferred embodiment such 
transforms result in a skewing or tilting of the image at angle 
4> and stretching of the image by a factor f, as given by the 
following general equation: 



x'-.t+y/-cos(<ji) 
and 



EQ. 1 



[0045] where (x, y) are the original image coordinates and 
(x', y') are the equivalent transformed or distorted image 

coordinates. 

[0046] More importantly, if the values for the parameters 
(t)and f can be determined then the image distortion effects 
can be corrected for by applying the reverse transformation 
of EQ. 1 as given by: 



EQ. 2 



[0047] where, (x, y) are the original image coordinates and 
(x\ y*) are the equivalent transformed or distorted image 
cooidinates. 

[0048] In turn, from geometric principles, the parameters 
(j)and f are determined by the orientation of the illuminating 
- x-ray beam relative to the surface of the radiation detector. 
Specifically, the skew angle, ((), is equal to the beam azi- 
muthal angle and the stretch factor, f , in the far field depends 
upon the inclination angle, 6, of the beam relative to a 
perpendicular to the surface of the radiation detector by the 
relationship: 



1 

' cos(fl) 



EQ. 3 



[0049] However, it can be readily seen that the angle, ^ 
and factor, f in EQ. 2 can be directly obtained from the 
extent, and direction of elongation of the image of the 

^spherical fiduciary element captured by the sensor, which is 
essentially the shadow of the fiduciary element cast upon the 
surfa<» of the radiation detector by the illuminating Xrray 

-beam. 
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[0050] FIGS. 3a and 3b arc idealized schematics illustrat- 
ing the geometric arrangement of a fiduciary element 114 
with respect to the surface of the radiation detector 118 and 
the illuminating x-ray beam and the resulting image. Note 
that if the direction of the x-ray beam is perpendicular to the 
surface of the surface of the radiation detector, the image 
will be a circle. FIG. 3b shows the relationship between the 
beam inclination angle, 6, and the length of the major axis, 
a, under the far field assumption. The relationship defined in 
EQ. 3 can readily be deduced from FIG. 3^. Note that the 
length of the minor axis, b, is the same as the diameter of the 
fiduciary element. 

[0051] From FIG. 3a it can be seen that if the beam is not 
perpendicular to the surface of the radiation detector then the 
shadow will be elongated and slanted or^skewed al the same 
angle as the beam azimuthal angle, ((>. Also, under the far 
field assumption^ the shape will be an ellipse, where the 
^length of the minor axis, b, will be the same as the diameter 
oMherfiduciary element and the length of the major axis, a, 
will be loiter than b by the factor, f, hence: 




EQ. 4 



[0052] . Evca in the case of a non-planar parallel x-ray 
Warn, more specifically, a divergent beam, the above deri- 
vation for f is robust. The beam divergence, here designated 
by Ae, for most intraoral dental applications is typically less 
than 5 .degrees. The error for f, defined as f(e*Aer4 and 
denoted as Af, can be calculated in terms of the beam 
dive^ence. More specifically, it can be shown from EQ. 3 
that the fractional error for f, defined as Af/f, is approxi-j 
mately: . ' ~ - 

r - cos(g) J . EQ. 5 



[0053] For S-M) degrees and A0=5 degrees the fractional 
error is less than 0.004 or 0.4 percent. For values of Oas large 
as 30 degrees and AG=5 degrees the fractional error is less 
then 0.06 or 6 percent. 

[0054] In addition, the error in the estimate of the skew or 
azimuth angle, <(), is not affected by beam divergence, 

[0055] In practice, the parameters a, b and (]), for the 
elliptical images of the fiduciary elements ellipse can be 
extracted from the image data using image processing 
techniques and numerical methods. In the preferred embodi- 
ment this analysis is performed on a digital representation of 
the radiograph. Note that the actual image-capture technol- 
ogy for a digitized radiograph is not the subject of the 
present invention. Such images may be obtained indirectly 
by scanning conventional film radiographs with currently 
available digital fihn scanners or obtained directly using 
existing digital radiographic image capmre technology such 
as storage phosphor plates or CCD/scinlillalor or CMOS/ 
scintillator combinations. 

[0056] The processing of image data in the preferred 
embodiment of the present invention is accomplished in two 



steps. First, the region of the digitized radiographic image 
that includes a fiduciary element is isolated based on its 
known position on the sensor. An example of such an image 
is shown in FIG. Sa for a spherical fiduciary element. Note 
the elliptical shape of the image due to the non-perpendicu- 
larity of the x-ray beam relative to the surface of the 
radiation detector capturing the image. 

[0057] Next, the edges of the ellipse are emphasized using 
gradient filtering (e.g. Sobel filter) as illustrated in FIG. Sb. 
An intensity threshold is applied to the filtered image to 
extract image points in the edge region of the ellipse, and 
then these points arc thinned or eroded using standard 
numerical techniques in order to more precisely define the 
edge region. These extracted points are shown in FIG. 5c. 
These points are used in the derivation of the curve fitting 
parameters. To this end, the set of extracted edge points 
(white points) in FIG. 5c are tabulated in terms of their pixel 
coordinates (i.e. x and y pixel coordinates). 

[0058] The coordinates of the tabulated points, are applied 
to the fitting of an equation for a generalized ellipse oriented 
at angle, <t», with major axis length, a, and minor axis length, 
b, and centered about (Xq, yo) as given by: 

r(jc-xo)-cos(<^)-l^ r(jc-JCo)-siii(^) + l^ ^Q- ^ 
I t>'->'o)-sin(^) 1 , 1 (y->-o)-cos(») J ' . 
P P = ^ 



[0059] Note that EQ. 6 is not the same as the image 
distortion transformation due to beam non-perpendicularity, 
but rather just the equation of a rotated ellipse used to 
estimate, the orientation parameters of interest, <j), a, b. 

[0060] The best estimate for the unknown parameters, <(), 
a, b, Xq, and yo are derived through least squares error 
minimization, minimizing the summed squared error 
between the coordinates of the tabulated extracted feature 
points and their corresponding equivalent points on the 
parameterized ellipse given by EQ. 6. In the preferred 
embodiment, an initial estimate is made for the unknown 
parameters based on an undistorted image and correspond- 
ing pairs of points are found by taking the nearest point on 
the parameterized ellipse to each extracted feature point. 
This approach, when applied iteratively, rapidly converges 
to an error minimum. In general, all such methods of initial 
parameter approximation and method of pairing correspond- 
ing points is claimed. 

[0061] Once the best estimate for the unknown param- 
eters, <|), a, b, Xq, and yo are derived one can then use the 
estimate of the parameters a and b to calculate the factor f 
through EQ. 4 

[0062] Note that in this specific example, the expression is 
not linear in all the parameters, <|), a, b, Xq, and yo, thus linear 
least squares fitting methods do not apply. However, numer- 
ous non-linear methods exist, such as iterative application of 
gradient descent error minimization. Also note, the use of 
other shapes for fiduciary elements, in other embodiments of 
the present invention may be chosen that offer more direct 
solutions to the orientation parameters. Finally, other meth- 
ods besides least squares error minimization may be applied 
to derive estimates of the beam orientation parameters. The 
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present invention is intended to cover all such fiduciary 
element configurations and methods of deriving orientation 
parameters from the image of such elements. 

[0063] The final step is to apply the corrective transfor- 
mation in EQ. 2 to the digitized image based on the derived 
values ()> and f. 

[0064] As an example of the effectiveness of the tech- 
nique, FIG. 5a shows an example of a radiographic image 
&om a bi-planar sensor with imbedded 1 mm spherical 
fiduciary elements. Note the distortion, specifically elonga- 
tion, of the roots of each maxillary tooth and the corre- 
sponding elongation of the fiduciary elements due to exces- 
sive x-ray beam inclination. 

[0065] FIG. Sb shows the result of correction of the image 
in FIG. 5fl using the extent of elongation of the fiduciary 
elements to correct the image. Note with the corrected image 
the fiduciary elements are restored to a circular shape. 

[0066] In view of the above, it will be seen that the several 
objects of the invention are achieved and other advantageous 
results attained- 

[0067] As various changes could be made in the above 
construction and methods without departing from the scope 
of the invention, it is intended that all matter contained in the 
above description or shown in the accompanying drawings 
shall be interpreted as illustrative and not in a limiting sense. 

We claim: 

1. A radiation sensor for use with a digital radiography 
imaging system for intraoral placement in a mouth of a 
patient for production of radiographs of teeth and supporting 
structures, said sensor comprising: 

a housing containing at least one generally planar, radia- 
tion detector providing an imaging surface oriented 
toward a radiation source, said housing conforming to 
the anatomic curvatures of the human maxillary and 
mandibular arches of the average patient; and, 

at least one radio-opaque fiduciary element of known 
shape, size and location positioned intermediate the 
radiation source and the surface of said at least one 
radiation detector, said at least one fiduciary element 
casting a projected image on said at least one radiation 
detector when illuminated by a radiation source. 

2. A radiation sensor as in claim 1 wherein said housing 
contains at least two, generaUy planar, radiation detectors 
abutting at a non-zero angle to form a faceted, generally 
contiguous imaging surface oriented toward the radiation 
source, said angle selected to conform said housing to the 
anatomic curvaUires of the human maxillary and mandibular 
arches of the average patient. 

3. A radiation sensor as in claim 2 wherein each adjoining 
pair of said generally planar detectors abut one another al a 
fixed angle. 

4. A radiation sensor as in claim 2 wherein each pair of 
adjoining generally planar radiation detectors are flexibly 
joined so thai the angle at which they abut can be changed 
to conform said housing to the anatomic curvatures of the 
human maxillary and mandibular arches of the average 
patient. 

5. A radiation sensor as in claim 2 wherein said fiduciary 
element is a sphere. 



6. A radiation sensor as in claim 1 wherein said fiduciary 
element has a shape whose projected image on the surfaces 
of said radiation detectors includes at least two intersecting •'^ 
line segments, " 

7. A radiation sensor as in claim 1 further comprising at 
least one radio-opaque fiduciary element of known shape, 
size and location embedded on, in or under said housing 
intermediate the radiation source and the surface of each of 
said at least one radiation detectors, said at least one 
fiduciary element casting a projected image on said at least 
one radiation detectors when illuminated by an x-ray source. 

8. A radiation sensor as in claim 1 wherein said housing 
further having a holding tab protruding therefrom for reten- 
tion ttelween the teeth for holding the radiation detectors in 
a fixed position in the patient's mouth. 

9. A radiation sensor for use with a digital radiography 
imaging system for intraoral placement in a mouth of a 
patient for production of radiographs of teeth and supporting 
structures, said sensor comprising: 

a housing containing at least two, generally planar, radia- 
tion detectors abutting al a non-zero angle to form a 
faceted, generally contiguous imaging surface oriented 
toward the radiation source, said angle selected to 
conform said housing to the anatomic curvatures of the 
human maxillary and mandibular arches of the average 
patient, said teeth and supporting structures casting a 
projected image on said at least two radiation detectors 
when illuminated by a radiation source. 

10. A radiation sensor as in claim 9 wherein each adjoin- 
ing pair of said generally planar detectors abut one another 
at a fixed angle. 

11. A radiation sensor as in claim 9 wherein each pair of 
adjoining generally planar radiation detectors are flexibly 
joined so that the angle al which they abut can be changed 
to conform said housing to the anatomic curvatures of the 
human maxillary and mandibular arches of the average 
patient. 

12. A radiation sensor as in claim 9 wherein said housing 
further having a holding tab protmding therefrom for reten- 
tion between the teeth for holding the radiation detectors in 
a fixed position in the patient's mouth. 

13. A method for correcting distortions in a radiation 
sensor used wi± a digital radiography imaging system for 
intraoral placement in a mouth of a patient for production of 
radiographs of teeth and supporting structures, said method 
comprising the steps of: 

Providing a housing containing at least one generaUy 
planar, radiation detector providing an imaging surface 
oriented toward a radiation source, said housing con- 
forming to the anatomic curvatures of the human 
maxillary and mandibular arches of the average patient; 

Placing at least one radio-opaque fiduciary element of 
known shape, size and location positioned intermediate 
a radiation source and the surface of said at least one 
radiation detector, said at least one fiduciary element 
casting a projected image on said al least one radiation 
detector when illuminated by the radiation source. 

Exposing said at least one radiation detector and said at 
least one fiduciary element to the radiation source to 
project an image of said at least one fiduciary element 
and said teeth and supporting structures onto the sur- 
face of said at least one radiation detector; 
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Capturing and digitizing the data representing the pro- 
jected image of the fiduciary element, teeth and sup- 
porting structures images on the surface of said at least 
one radiation detector produced by the radiation detec- 
fvtor. 



"Analyzing the digitized image data to determine" the j 
distortion of the projected fiduciary shape onto the 
surface of the at least one radiation detectors due to the 
non-perpendicularity of the radiation source with 
^respect to the surface of the at least one radiation 
' detectors from that of a ideal fiduciary-image projected 
onto the surfaces of the at least one radiation detectors J 
: by exposure to a radiation source perpendicular to the 
^surface of the at least one radiation detectors, and 
determining a corrective transformation that transforms 
the distorted projected fiduciary shape to that of the 
ideal fiduciary shape; and, 

Appfying the corrective transfonnation to theTeinaining 
' digitized image data in order to transform the distorted 
projected image of the teeth and supporting struciuies 
to Ih at ~oi an ideal projected image of the teeth and 
supporting structures. — 



14. A method as in claim 13 further including the step of: 

Placing a radio-opaque fiduciary element of known shape, 
size and location intermediate the radiation source and 
the surface of each one of said at least one radiation 
detector 

15. A method as in claim 14 wherein the step of placing 
at least one radio-opaque fiduciary element of known shape, 
size and location intermediate the radiation source and the 
surface of each one of said at least one radiation detectors, 
further includes the step of: 

Embedding on, or placing in or under said housing at least 
; one radio-opaque fiduciary element of known shape, 
size and location intermediate the radiation source and 
the surfaces of each of said at least one radiation 
detectors, said at least one fiduciary element casting a 
/ projected image on the surfaces of said at least one 
radiation detectors. 
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A method to calculate angular disparities between object and 
receptor in ^paralleling technique' 

RKW Schulze* ' and B d'Hoedt' 

^Department of Oral Surgery , Dental School. Johannes Gutenberg University, Mainz, Germany 



Objective: To develop an image analysis method for calculation of angular disparities between 
an object, temporarily equipped with a reference system, and a radiographic receptor. 
Materials and methods: A mathematical method based on a reference system containing two 
metallic spheres is developed, allowing calculation of inclination between the inter-spherical 
axis and the digital image receptor using image features. Experimental evaluation was done in 
standardized projection geometry for two sphere sizes at four randomly chosen inclinations per 
size, with each radiograph assessed three times. Truth was assessed threefold from photographs 
obtained at each inclination. 

Results: Mean standard deviation between single assessments was 2.6°. Significant differences 
{P Maloney/Rastogy = 0.00) were found between absolute values of truth and calculated values 
(mean: —0.9°; range: —6.0°; 3.6°), indicating a significant lack of accuracy. 
Conclusions: Although so far not sufficiently accurate, the method yields information relevant 
for correction of distortion in intra-oral radiology. 

Dentomaxillofacial Radiology (2002) 31, 32-38. DOI: 10.1038/sj/dmfr/4600657 

Keywords: image processing, computer-assisted; computer simulation; radiography, dental 



Introduction 

Geometric discrepancy between serial X-rays may 
indicate differences that do not reflect anatomical 
changes.' Standardization of the imaging geometry is 
therefore a prerequisite of utmost importance, if 
accurate measurement of such changes is required. 
According to a fundamental study of van Aken the 
most desirable projection is obtained, when (1) parallel 
rays are used due to infinite distance between focal 
spot and object, and (2) the X-rays meet both film and 
object at right angles.^ Since criterion (1) cannot be 
achieved as the focus is always at a finite distance to 
the object, the paralleling technique was introduced in 
dental radiology as a good compromise. While a rigid 
holding device provides a perpendicular alignment of 
the central ray relative to the receptor, true parallelism 
between receptor and object is hardly feasible in a 
clinical environment, particularly in the upper jaw. 
Often the object and receptor plane are not parallel, 
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resulting in a non-affine and irreversible distortion of 
the image as a function of degree of actual 
inclination.^'* Angular differences between two expo- 
sures can be derived mathematically from images, if 
radiographically detectable landmarks of known 
dimensions are present and the geometric relation 
between focus and object is fixed. However, no 
method exists to calculate the actual angulation of the 
receptor plane relative to the object in one radio- 
graphic image. In tuned aperture computed tomogra- 
phy (TACT) metallic spheres can be used as fiducial 
markers to obtain the center of gravity for reconstruc- 
tion of the TACT image.* In addition, six or less 
spherical references can be used to determine the 
geometrical positions of the X-ray source, if the object 
is equipped with references and a back projection 
algorithm used.^ In contrast to this method, our 
approach is aiming to determine angular disparities in 
paralleling X-ray projection by quantitative analysis of 
the elliptical shape of two spherical references and by 
location of their shadows. 

We based our investigation on the fact that 
magnification is always directly proportional to the 
distance of an object from the projection plane."* 
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Hence, by means of a reference object of known 
dimensions, the position of the object within an 
otherwise fixed projection system can be derived from 
its magnification. We believe that the angular relation- 
ship of object to receptor may be obtained from 
quantitative image analysis using two reference spheres 
temporarily attached to an object with their main axis 
aligned with the object. 

The aim of this investigation was to calculate the 
angular disparity of an object relative to the receptor, 
based on the application of reference objects and a 
perpendicular orientation of the central X-ray to the 
image receptor. 



Materials and methods 

Mathematical theory 

The following method is restricted only to the 
paralleling technique and is based on the fact that 
transmission radiography of an object with rotational 
symmetry, such as a sphere, will only be distorted in 
one direction, depending on its actual position in the 
projection system. If the central ray passes through the 
center of the sphere the resulting image will be a disc, 
but deviation from this path will produce an elliptical 
image of the sphere. The magnitude of the distortion is 
a function of three factors: (1) the distance between the 
focal spo t and the sphere as measured along the central 
ray {FX,)\ (2) the distance between the sphere and 
receptor as measured along the central ray {XyO')\ and 
(3) the shortest distance of the sphere from the central 
ray of the beam (d) (Figure 1). If the dimensio ns o f the 
sphere and the focus to receptor distance {FO') are 
known, its actual position in a 3-D coordinate system 
with the X-axis aligned to the central ray can therefore 
be derived directly from the image. When a sec ond 
sphere is applied, the inter-sphere axis {M1M2) as 
defined by the line connecting both centers determines 




>ibcal8}iot(F) 



Figure I Underlying projection geometry resulting in one-dimen- 
sional distortion of objects with rotational symmetry (spheres) in 
transmission radiography. Note the elliptical distortion caused if the 
centr al X-ray is not passing thro ugh the center of the sphere. 

= object to receptor distance; FA'„ = focus to object distance; 
d = shortest distance between center of sphere and central X-ray 



the actual inclination relative to a receptor axis of 
interest. 

Our calculations are based on the following 
assumptions: 

(1) the focal spot (F: Xp, 0, 0) is approximated by a 
point, located on the x-axis of the 3-D coordinate 
system (X, Y, Z) with the central ray aligned with 
the AT-axis meeting the receptor surface at point 
zero (O': 0; 0; 0), 

(2) X-rays are described by linear rays, 

(3) y- and z-axes define the horizontal and vertical 
longitudinal axes of the flat receptor, 

(4) the inclination is calculated in relation to an axis of 
interes t of the receptor (here: vertical axis (z)), with 
M1M2 being orientated relatively parallel to this 
axis in the >^-z-plane. 

(5) the position of the central ray on the rece ptor 
surface and the focus to receptor distance (FO^) are 
known. 

(6) the position of the focal spot is not coHnear with 
the centers of the reference spheres. 

All subjects of interest are graphically illustrated in 
Figure 2. For simplification reasons, the following 
explanations are presented in a generalized manner 
with the index i replacing the number of the sphere (1, 
2), The upper sphere is referred to as sphere No. 1, the 
lower sphere as sphere No. 2. Let the objects located in 
the image plane be indicated by a superscript dot-dash 
('). The centers of the spheres are denoted by Mi, their 
images as M'j, respectively. R'j and S'j represent the outer 
and inner intersections of the imag e ellipse with the line 
passing through O' and Mj' (O'A/J). The angle ol, is the 
outer angle defined by R'i, F and S';, the inner angle Pi by 
S'i, F and O'. We obtain the angles a-, and fi, from: 



a; = -==- — Pi = arctan 
O'F ' ' 



arctan =1= (1) 
O'F O'F ^ * 



We now need to calculate the distance A/,F of the 
center of the sphere (Mi) from the focal spot F. By 
introduction of an auxilia ry plane defined by O'M'i and 
F we can calculate MiF as follows: 

sm(f) 

where x, is the radius of the sphere i. Since we are 
interested in the distance of Mii relative to the vertical 
axis (z) of the receptor, we have to find the projection 
Vi (Xvi, Zvi) of Mi onto a plane defined by .v- and z-axis 
of the coordinate system. To obtain these coordinates, 
we first have to construct the image of Vi (V'j) by 
drawing a perpendicular line to the vertical axis z 
through M'i. We then need to calculate th e ang le e 
between the X-ray passing through Mi and MiXy. in 
terms of: 



(3) 
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Figure 2 Geometrical construction for mathematical approach. The point of incidence of the central X-ray on the receptor (point zero: 
O'(0,0,0)) has to be known. R and S denote the outer and inner interception between the ellipse of the sphere image and the line through O' and 
the center of the elliptical image of each sphere. Detailed description of all other distances and landmarks is given in the text 



Xvi is derived from: 
with: 



(4) 



FXv, = A/,Fsin e (5) 

To clarify the underlying mathematics of the following 
steps a plane projection extracted from Figure 2 is 
presented in Figure 3. Zvi is obtained from the normal 
equation (eq.) of linear functions in terms of: 



VO' 



(6) 



Please note, that the vertical position of V, and V'2 to 
point zero (O') is essential for further calcu lation s, 
since it determines the direction of inclination. V\0' is 
signe d positive, if V is located above and, vice versa, 
V\0' is signed negative, if V' is located below point 
zero. Since we now have defined the X and Z 
coordinates (Xvi, ZvO of the projection of each 
spheres' actual position in the projection system onto 
the vertical receptor pl ane, we can now calculate the 
actual inclination y of M1M2 relative to the vertical 
axis of the receptor from: 



y = arctan 



Zv\ — Zv2 



(7) 



with y signed positive, if Xvi^Xv2, meaning that 
MxMi is inclined with its upper part towards the 
focal plane. Vice versa y will be signed negative, if 




1 



Figure 3 Two-dimensional projection extracted from Figure 2 for 
clarification of equations (6) and (7), which are necessary to obtain 
the actual angulation relative to the receptor. V,(Xvi; Zvi)/V2 (Xv2; 
Zv2) represent the projection of the centers of the reference spheres 
(M1/M2) on the vertical plane defined by focal spot (F) and vertical 
receptor axis. For further descriptions see Figures I and 2 



Xvi<Xv2, indicating an angulation of A/1A/2 with its 
top towards the receptor plane. 



Experimental set-up 

Experimental evaluation of the mathematical algorithm 
was performed by means of a projection system 
fulfilling the criteria as described above (Figure 4). 
Two steel balls (diameter: 3.15 mm) were embedded 
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13 mm apart from each other in clear autopolymeriz- 
ing resin in a small clear acrylic tube (length = 30 mm; 
diameter = 9.0 mm), with a minimum of 5 mm resin 
covering the balls on both sides. The same procedure 
was carried out with two additional steel balls of 
5.00 mm diameter. Both tubes were then cut at one 
end in such a way, that an oblique bevel resulted 
which was then glued to the top of a screw (M8; 
length = 20 mm). A three-sided open box 
(length = 14 cm; width=12cm; height=12cm) made 
of wood was constructed with walls standing 
perpendicularly to each other. In the following the 
long vertical wall (14 cm) will be referred to as back 
wall, and the smaller (12 cm) vertical wall as front 
wall. On the ground surface of the box a plastic 
positioning ring of a commercially available film 
holding device (RWT, model anterior, Kentzler- 
Kaschner Dental GmbH, Germany) was fixed with 
the ring being in parallel position to the front wall at 
a distance of 48.0 mm. This positioning ring has a flat 
front side and a back side with a rectangular notch 
exactly fitting to the end-piece of the radiographic 
tube (Hehodent MD, Sirona Dental Systems GmbH, 




Figure 4 Photo of experimental set-up of the situation encoded as 
image 2A. The left arrow is indicating the CCD-sensor, the right 
arrow the tube containing the two reference spheres presenting a 
forward (negative) angulation relative to the CCD-receptor. Light 
refraction of the clear resin within the cylindrical tube functioning as 
convex lens causes horizontal widening of the spheres in the 
photography. All relevant distances and dimensions are listed in 
Table I 



Bensheim, Germany). A disc of exactly the size of the 
outer diameter of the positioning ring (74.0 mm) made 
of thin (0.1 mm) clear plastic foil and containing a 
small metal ball (marking sphere; diameter = 0.80 mm) 
fixed in its center was concentrically glued to the front 
side of the ring. At a distance of 21.0 mm from the 
front wall the corresponding nut of the screw was 
mounted on the ground surface in such a way, that it 
was positioned slightly laterally to the central X-ray 
when the tube was connected to the positioning ring. 
By turning the screw in the corresponding nut various 
inclinations of the two-ball-system relative to the 
receptor could be produced. To obtain information 
on the actual inclination of the reference device 
relative to the vertical axis of the receptor, milli- 
meter-scaled paper was glued on to the back wall of 
the box aligned with the vertical axis of the front wall. 
A CCD-image receptor (40 x 24 x 6 mm) providing 
high resolution (Trophy RVG 5, Trophy Radiologie 
S.A., Marne la Valle, France) was rigidly mounted on 
to the front wall with its active surface in beam 
direction and the vertical ( = long) edge aligned 
parallel to the back wall of the box. Containing an 
active array of 30 x 20 mm this sensor yields a 
resolution of 1 572 864 pixels with each pixel being 
19.5 X 19.5 /im^ One option included in the proprie- 
tary software of the manufacturer (Trophy Windows 
Version 4.1k) is 1:1 pixel-binning, meaning that one 
sensor pixel will be represented by one image pixel. 
The sensor was interfaced to a personal computer 
(Intel Pentium III-MMX). After rigid fixation of the 
tube in the corresponding notches of the positioning 
ring by means of adhesive tape a series of four 
radiographs with each tube positioned in its nut was 
taken (60 kVp, 7 mA, 0,10 s) with different randomly 
chosen inclinations of the two-sphere system. For each 
inclination a photograph was obtained with a digital 
camera (FinePix SI Pro, Fuji Photo Film (Europe) 
GmbH, Dusseldorf, Germany) positioned approxi- 
mately I m distant to the box so that its optical axis 
was aligned with a theoretical line parallel to the front 
wall of the box through the center of the screw. Due 
to this set-up minimal distortion of the photos (jpcg- 
files) was guaranteed. Radiograph (tiff-file) and 
corresponding photography (jpeg-file) were exported 
and each pair was encoded for subsequent image 
evaluation. All relevant projection parameters are 
listed in Table 1. 



Table 1 Projection and exposure parameters for the experimental evaluation 



Focus to receptor distance {FO') 261.15 mm 

Distance between marking sphere and receptor 45.30 mm 

Diameter of the two reference spheres (ri, rj) rj: 5.00 mm, 4 exposures 

r2: 3.15 mm, 4 exposures 
Distance between reference spheres ri:l3.35 mm 

r2:l6.35 mm 

Pixel size of receptor RVG 5 0 0195 mm^ 

Exposure parameters 60 kVp, 7 mA, O.IQ s 
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Image evaluation 

Radiographs were examined by image-editing software 
(Adobe Photoshop 5.0, Adobe Software Inc. Mountain 
View, CA, USA) and then evaluated on a 1:1 display 
monitor (1028 x728 pixels) with each image pixel 
represented by one monitor pixel (Figure 5). The 
actual length of an image path can be calculated from 
the number of pixels and their length (19.5 /im) by the 
Pythagorean theorem. After brightness and contrast 
enhancement of the radiograph, the center of the 
circular image of the marking sphere (O': 0; 0; 0) was 
constructed by bisecting the horizontal and vertical 
maximum diameter. Subsequently, tangents to both 
sides of the images of each reference sphere were drawn 
and the resulting distance between the tangent points 
at each sphere bisected. By connecting O' to this 




Figure 5 Geometrical construction on a radiographic image 
obtained with the experimental set-up, required to extract all 
relevant information on distortion for calculation of actual 
angulation. For better understanding landmark points as introduced 
in Figure 2 are included in the image. The arrow is indicating the 
circular image of marking sphere, with its center representing the 
point of incidence of central X-ray (point zero: O' (0; 0; 0)) 



bisecting point we obtain the axis along which the 
sphere is distorted, and consequently, the construction 
points Ri' and S/. By again bisecting the distance 
between and S/ we obtain point M{, The 

intersection between a line through point O' 
parallel to the long image axis and a perpendicular 
dropped to this line through Mi' defines Vj' (Figure 5). 

Now the dis tances requi red for the calculation 

(Rx'O', Ri'O', Sx'0\ S2'0', V\0\ V'^O') were mea- 
sured to the nearest pixel and the values transferred in 
a spreadsheet software (Microsoft Excel 2000), where 
the resulting lengths were calculated. Again these 
values were transferred into another spreadsheet 
incorporating the described algorithm, and y was 
calculated for each radiograph. The whole construc- 
tion and measurement procedure was repeated three- 
fold by one examiner (RS), with the mean value 
representing ^calculated inclination'. 

The corresponding photographs were also evaluated 
with Photoshop software in 1 : 1 display. By counting 
horizontal and vertical distance between the spheres in 
millimeters on the millime ter-scaled paper the actual 
'true' inclination of M\M2 relative to the receptor was 
calculated by the equation: 



y = arc tan 



horizontal distance [mm) 
vertical distance [mm] 



This procedure was also repeated threefold by one 
observer (RS) and the resulting mean is referred to as 
*true inclination'. 



Results 

Considerable standard deviations (s.d.) within the three 
replicate assessments were found, with a mean s.d. of 




Figure 6 Modified Bland- Altman-ptot of observed (true) and 
calculated values. Differences between truth and calculated angula- 
tions (y-axis) are plotted against ground truth on the ;c-axis, with 
zero difference indicating maximum accuracy. Limits of tolerance are 
included in terms of mean difference (dark horizontal line) ±2 
standard deviations (gray shaded horizontal lines). The plot indicates 
limited accuracy of the method with more over- than under- 
estimations of the calculated values 
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2.6°. True' and ^calculated' values are compared in 
Table 2. To obtain information on the actual direction 
of bias for the calculated differences between truth and 
radiographically derived values, this calculation is based 
on absolute values. Otherwise we would have calculated 
only the magnitude of error, but not its direction with 
respect to over- or underestimation. Although on 
average calculated values differed from truth only by 
— 0.9° (s.d. = 3.5°), this was highly significant 
(P = 0.000) as indicated by the Maloney/Rastogi test." 
The root mean square of the actual differences was 3.4°, 
also indicating the lack of accuracy. 

Differences between truth and calculated values are 
graphically presented as modified Bland-Altman plots 
(Figure 6). Differences between individual pairs are 
plotted against truth, with zero indicating maximum 
accuracy. Limits of tolerance (gray-shaded horizontal 
lines) are included in terms of mean difference (dark 
horizontal line) ±2 s.d. No value exceeds these limits, 
however, considerable deviations of single values from 
truth (range: —6.0° -3.6°) are obvious from the plot, 
mainly indicating overestimation of truth. 

Based on the differences between absolute values, no 
significant influence of sphere size on accuracy was 
computed (Pwncoxon = 0.7l5). 



Discussion 

Our investigation deals with calculation of actual 
degree of angular disparity between object and image 
receptor from intra-oral radiographs acquired in 
*paralleling-technique'. Such disparities are common, 
since the angle is often difficult to control in a clinical 
situation. Considerable distortion of the images under 
examination may result from angular discrepancies.-^ 

Our method uses two radio-dense objects with 
rotational sym metry (metallic spheres), with the inter- 
spherica! axis (AfiA/2) as defined by the connecting line 
between the centers of the spheres being the axis under 
examination. Prerequisites for the applicability of the 
technique are true perpendicular orientation of the 
central X-ray relative to the receptor in combination 



Table 2 Photographically derived 'true' and radiographically 
computed values for all images of both series. Note, that for 
calculation of the difference between truth and calculated value the 
signs (directions) of the angles were omitted to obtain information on 
the direction of the bias 



Series/ image 
No. 


True Calculated 
inclination O inclination O ± s.d. 


Truth - calculated 
value 0 


IX/2A (r,) 


-16.6 


-18.5±1.8 


-1.9 


1X/2B (r.) 


-7.8 


-13.8±2.5 


-6.0 


IX/2C (r.) 


19.5 


15.9±l.3 


3.6 


IX/2D (r.) 


-6.5 


-3.2 ±3.7 


3.3 


X/2A (r^) 


-2.2 


-2.8±3.7 


-0.6 


X/2B (rj) 


16.6 


I5.1±5.5 


1.5 


X/2C (r2) 


-5.4 


-9.4 ±0.9 


-4.0 


X/2D (r2) 


-19.0 


-22.4 ±1.5 


-3.4 



s.d., standard deviation; r, radius of reference spheres (rp 5.00 mm; 
3.15 mm) 



with a known point of incidence (0, 0, 0) and a known 
focus-receptor distance. In principal, the method is based 
on measurement of magnification of the spheres, since 
their shadows are elliptically distorted only in one 
dimension due to their symmetrical design. From these 
data the position of the respective sphere within the 
projection system is calculated, and the positions of both 
spheres relative to a specified receptor axis determine the 
actual inclination. Although any desired axis within the 
receptor plane could apply for the investigation, we 
selected the axis parallel to the long object axis, since 
angular discrepancies will most probably occur in this 
direction. If the two spheres are temporarily fixed in a 
known position to the object of interest (tooth, dental 
implant), the actual degree and direction of inclination 
relative to the receptor plane could be calculated. Our 
method could be used for postoperative monitoring of 
peri-implant bone height. It is a well known fact, that 
exact measurements require strict parallelism between 
fixture and receptor axis.'^ '^ This prerequisite could, 
however, be of less importance if information on the 
actual inclination was derived from the radiograph, 
making a posteriori correction possible. If for instance 
the reference spheres were implemented in a temporary 
abutment fixed in axial direction to the implant during 
radiographic exposure, this information could be 
obtained. A distorted linear distance within such an 
object, such as the vertical long edge of an implant, could 
then be corrected for distortion. For pre-implant 
planning, the spheres could also be placed into pre- 
surgical stents with M1M2 representing the desired 
implant direction. From the calculated angulation linear 
bony dimensions in axial direction of the implant could be 
more accurately corrected for distortion than with any 
other reference (single sphere, grid). In periodontal 
assessments of marginal bone height, the spheres could 
be attached to the proximal surface of a single tooth 
aligned with its main axis. Linear bone height measure- 
ments in this direction could then also be corrected 
analogously. In addition, information provided by the 
method could be used in subtraction radiology to limit the 
degrees of freedom between subsequent radiographs by 
definition of the position and angulation of an object.'^ 
Precise calculation, however, requires information of 
angulation of the object plane relative to the receptor 
plane in horizontal and vertical direction, which could be 
obtained from the method by means of a triangle 
containing three spheres placed in the object plane. We 
are currently working on the mathematical algorithm for 
correction of linear distortion from such angulations. 

Our experimental data show, that, in principal it is 
possible to derive the actual degree of angular disparity 
from the two-sphere reference system. However, with the 
applied set-up both accuracy as well as reproducibility 
are clearly not sufficient for clinical use. The trend of 
overestimation, however, should be seen in the context of 
the small sample size (eight exposures) and should 
therefore not be considered as methodological bias. 
Obviously, our assessments show low reproducibility as 
indicated by a mean standard deviation of 2.6° (mean) 
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between single assessments. The design of the experi- 
mental set-up is one fundamental origin of error. 
Presumably, marking of the essentially important point 
of incidence of the central ray (point zero) with the steel- 
sphere mounted in the center of the tube aperture was 
not sufficiently accurate, and in addition, its shadow was 
blurred due to the long sphere-receptor distance 
(45.30 mm, Table I). Technically this problem could be 
resolved by means of a rigid holding device with a mould 
for the digital receptor centered relative to the central X- 
ray in such a way, that the central receptor pixel 
represents point zero (0,0,0). The device must also 
ensure a true perpendicular alignment of the central 
ray relative to the receptor and a normal cone length 
should be used. The latter is recommended since a long 
cone would result in less magnification and, conse- 
quently, in increasing method error. 

Another source of error is the assessment procedure of 
magnification and shadow location, which was performed 
by simple linear image measurements. Accuracy of the 
method is strongly depe ndin g on measurements of the 
dimensional changes with an approximate 

magnitude ranging between 0.01 mm and 0.1 mm, 
dependent on the actual position of the sphere within 
the system and on the diameter of the spheres. Digital 
images are limited with respect to the spatial resolution 
due to their discrete sampling character. Hence we needed 
a receptor with a small pixel size and the option to display 
each pixel separately (1:1 mode). The RVG 5 CCD- 
sensor provides a pixel size of 0.0195 mm^ and the 
software incorporates the option to display each pixel 
without pixel-binning. Consequently, changes exceeding 
0.0195 mm can theoretically be detected with this 
receptor, yielding sufficient resolution for our method. 
In theory, one would expect better accuracy for reference 
spheres with larger diameters, since absolute dimensional 
changes are larger with increasing diameter. For the small 
sample size, however, no such trend was evident from the 
data. At the same time, correct detection of those pixels 
that determine the end-points of the assessed distance is of 
fundamental importance for the accuracy of the 
procedure. We encountered two major problems during 
this procedure. First, the difficulty to select the correct 
pixel and second, to mark this pixel with the mouse-driven 
cross-shaped cursor. In addition, on a discrete pixel 
matrix, distances can only be defined by a multiple of 
whole pixels, with the result possibly being an odd 
number. Bisecting of such distances then produces *hair 
pixels, which then have to be assigned either to one or the 
other side of the distance assessed by subjective decision 
of the reader. Automated image enhancement and 
registration would facilitate standardization of the 
assessment problem. Since we have an a posteriori 
situation when calculating parameters based on informa- 
tion derived from a given radiograph, image registration 
features as described by Lehmann et aL apply, with the 
images of the spheres representing extrinsic features.'^ 
Automated detection of local contrast (defined as 
brightness difference between two neighbouring pixels) 
allows for accurate finding of simple structures as 



presented by the sphere images. By means of grayscale 
spreading and a predefined threshold those discs could be 
accurately detected, followed by an automated assess- 
ment of all required distances. In addition, ass essme nt of 
the area of the shado w ins tead of the distance R\S\ itself 
and then calculating R\S] from the principal axis of the 
ellipse could presumably enhance accuracy. These 
procedures could be implemented in any radiographic 
imaging software. We are currently initiating a software 
development for these purposes in our department. 

In conclusion, our data show that it is possible to 
derive information a posteriori on the actual inclination 
between a receptor axis of interest and an object 
temporarily equipped with two spheres. Refinements to 
the method may increase the accuracy to a clinically 
acceptable level and further research is required. 
Currently, no other method exists providing equal 
information on this geometrical parameter inducing 
distortion on radiographic images. 
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Determination of projection geometry from quantitative assessment of the 
distortion of spherical references in single-view projection radiography 
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A method is introduced, inferring the three-dimensional (3-D) location from the 2-D radiographic 
shadow of an opaque spherical reference body of known radius by considering its elliptical distor- 
tion, the 2-D shadow location and a known source-to-receptor distance. Three noncolHnear spheres 
fixed to a rigid object constitute all possible degrees of free€k>m. I.e., the entire 3-D imaging 
geometry. The method may be used (a) to detemaine the 3-D imaging geometry from a single 2-D 
view and (b) to correct for foreshortening of object distances coplanar with the plane defined by the 
sphere triplet. Apart from the mathematical background the article describes a smaJl feasibility 
experiment, performed with four different sphere diameters and a commercial dental ccd-receplor 
system (pixel length: 0.0195 mm). The mouse-cursor based image evaluation revealed an average 
underestimation of the critical depth- U-) coordinate decreasing with increasing radius (-30.3% for 
r=0,5 mm to 2.8% for r«2.5 mm). Intraobserver reliability (the standard deviation between three 
single cursor-based assessments) ranged between 0% and 8% of the actual true depth. The main 
source of the input error is associated with the assessment of the amount of elliptical distortion, 
where subpixel accuracy is demanded. Consequently, software-based automated image' evaluation is 
required using available methods for pattern recognition and point-spread correction. Provided 
sufficient accuracy, the method provides an important tool for foreshortening correction, depth 
assessment, motion analysis, and 3-D reconstruction from two or more 2-D views. O 2004 
American Association of Physicists in Medicine. [DOl: 10.1118/1.1796951] 



I. INTRODUCTION 

The spatial relation between x-ray sotirce, object, and imag- 
ing plane relative to one another during exposition, i.e., the 
projection or imaging geometry, fundamentally determines 
the two-dimensional (2-D) radiographic image formation. 
Since three-dimensional (3-0) reconstruction is commonly 
based on the determination of the imaging geometry, " its 
knowledge enables straightforward 3-D' reconstruction from 
at least two 2-D images obtained under different 
geometries.^*^ Knowing the imaging geometry is also helpful 
to calculate the magniftcation and foreshortening of objects 
appearing in the image. Such effects are very common, par- 
ticularly in dental intraoral radiography, where anatomical 
barriers (e.g., the hard palate) preclude a parallel positioning 
of object and image plane. Commonly, spherical reference 
markers are used for calibration,^"^^ preferably in combina- 
tion with a holding device'^'^ coupling the receptor to the 
source in a c-arm fashion. In medical radiography, metallic 
reference markers*****' or other objects of known dimensions, 
e.g., angiographic catheters,^** are used as calibration 
objects. Despite their frequent use in dental^"' ^ as well 
as medicar^'^' radiography or for tomosynthetic 
reconstruction,**' so far only magnification^* ^"^^ or the 
2-D shadow position**' of spherical reference bodies was 
considered, thereby neglecting the depth information inher- 
ent in the elliptical distortion. Depth information, however, is 
essential for the calculation of the imaging geometry. Hoff- 
mann et ai. introduced a method to determine the 3-D posi- 
tion of a general , sparse object from a single view by using 



a priori knowledge of relative marker positions. Reference 
spheres have already been used for 2-D assessment of the 
angular relationship between an object and a given receptor 
axis;^^ We aim at the extension of this method to the third 
dimension to obtain their 3-D position from a single projec- 
tion. It will shown that from quantitative exploitation of the 
2-D shadows of three noncollinear reference spheres fixed to 
an object under study, the imaging geometry of the exposure 
can be determined, if the source-to-receptor distance is 
Imown. A particular application of this method for the fore- 
shortening and magnification correction of objects coplanar 
with the reference plane determined by the sphere triplet will 
also be addressed. In addition, results of a small feasibility 
experiment are presented. 

II. METHODS 

A. Mathematical model 

7. Calculation of the 3-D reference sphere 

coordinates 

Consider a cone-beam projection in a Cartesian coordi- 
nate system (Fig. I), where the central x ray (aligned with 
the * axis), i.e., the orthogonal shortest distance between a 
spot source F(x/-,0,0) and a fiat y-z-receptor plane incides 
the latter in the origin O' (0,0,0). Let all 2-D images points 
be denoted by an apostrophe. The points of incidence of the 
diverging beam may be considered as an infinite number of 
concentric circles with center point O' (0,0,0). Thus, the 
maximum magnification on the 2-D image plane will occur 
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Fro. 1. Exposed from a focal spot 7^x^,0. 0), the maxi- 
mum tnagnificciion defining the major ellipse axis R'S* 
occurs along a Lme through the point of incidence 
O' (0,0,0) of the centnl x ray on the flat >r~2-image 
p\tMt..R*iO,ygt,2jt') ond S'iO^ys'.xs') represent the 
2-D images of the 3-D sphere tangent points 
R(xjt,yjt,ZK) and 5{x^,>'5,zj). They subtend an angle or 
betweai R' F S' denotes the angle S' F 0\ Note, 
that the image A/*(0,y^,x.v') of the desired 3-D sphere 
center point M{xj^,yij^Zj^) 6oe.% not coincide with the 
ellips e center, i^e^.does not bisect the long ellipse axis 



along vectors pointing frorn O- toward the periphery. We aim 
for the localization of a given sphere (radius: r), i e., its 
center of mass point JW(x^rf,yA^,2Af) frorii its shadow that will 
be an ellipse (for proof sec Ref.. 24), excjept for an orthogonal 
projection, where it will form a circular di^k. Thb projection 
lines tangent to the 3-D sphere tangent po ints R {xRiyR,Zji^ 
and 5(x5,3^5,z^) cast the major ellipse ajds i^ 
sion of which will riecessarUy intersect with due to the 
direction of 2-D miagnification. Thc.angle a (Fig. 1) between 
R*F and S* is given by 



tan p, 



where 



y3 « tan' 



(I) 



(2) 



Please note that the image AT .of the sphere's cehtcr jpoint Af 
does not coincide with the ccn^r of the- ellipse cast by the 
sphere. The distance between M and the . spot soiiixie is 
obtained from 



(3) 



From the theorem of intersecting lines we obtain the desired 
3-D coordinates of M from 



(4) 



(5) 



z^,|fiWlcos^| + ;9^ 



(6) 



The 3-D coordinates of the center points f^i* of thie 
three noncollinear spheres temporarily to a rigid object under 
study provides the 3-1) locations of three "object'* points. 
0'(6,0V0) may simply be, reconstructed from intersecting the 
lines extended through the long elliptical axes, or, in the 
presence of error, from the point ofclqsest approach between: 
tliese imes. Together with the Imbwn source-to-receptor dis-- 
taniae xj^* this information constitutes the spatial relation be^ 
tween the soiirc leaving the object three 

treii^l^dbnal and three rotati6ii|J degrees of freedonii Theses 
are sufficiently constitutejd* by the three points in, space; 
obtained fi^bm Eqs:^ 



2. Correction of foreshortening for 
withm-feferehcerplane ob^ 

If the noncollinear sphere triplet is placed coplanar .with 
an object G or a distance under examination (Fig. 2); the 
plane P defined by A^i, A/2» 



with obtained from the deteirninant 



x-x^ y-yx z-zi 
X2'Xi yi-yi Z2-Z1 
x^-xi yi-yi z^-zi 



= 0, 



(7) 



(8) 



may be rotated about the y axis at an arbitrar>' angle y and 
about the z. axis an angle e, respectively. By backprojecting a 
given point T {xj-^ ,yr »^r) ^ Q' onto P, we find its actual 
3-D position Ttxr^^'r.^r) e C from intersecting the line 
through FT with P, The equation of a line-plane 
intersection^* yields (for a derivation sec the Appendix): 
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Fio. 2. If the sphcre^triplet.is placed, coplanar with an object Q under ex- 
amination, any object- image point T{xr ^Vv .^r) ^ 6' pan be back- 
projected onto the reference plane /*, yielding, its tnic 3-0 position 



lining coordinatt 



The remaining coordinates arc given by 
and 



(9) 



(10) 



(U) 



B. Experimental eyaluatlon ' 

Four acrylic phantoms (thickness: 2.0 mm), each contain- 
ing three isometric steel spheres (0,5, 1 .0, 1 .5, and 2.5 mm) 
and two pieces of steeWire "objects" (diameter. 0;6 ram; 
lengths: 8.90 or 16.20 mm) with, needle-shaped , ends, were 
exposed separately on an optical, bench with, a fabrication 
accuracy of 0.1 to 0.2 mm. The point of incidence Q' (0 , 0 i O) 
of the central x ray coincided with the center of a dental: ccd 
receptor (Full Size, Sirona /Dental Systems GnabH; 3fen- 
sheim, Germany; physical pixel size: 0.0195 min 
X0.0195 mm) to within ±10 pixels. This design eliminated 
additional errors to be expected from distortion-based origin 
reconstruction from only minimally distorted shadows on the 
small-size (26 nun X 34 mm) sensor area. Since no image 
intensifier or fiber-optic taper was used and the sensor size 
was very small, negligible pincushion distortion is assumed. 
Each phantom was exposed at three time settings (0.06, 0:08, 
and 0.10 s) and eight different angulations (range: 
-57.7^-46.7'*) relative to the image receptor, occurring 
about the y axis and the z axis (Fig. 2). Within the series of 
0.08 s exposure time, scatter equivalents of one, three, or 
five wax plates (thickness: 1.5 mm) were employed in a step- 
wise fashion. Due to the construction of the optical bench, Xf 
ranged between 254.3 and 294,1 mm. For the mouse-driven 




Fig, 3. Using the mouse-driven cursor, tangents were manually drawn from 
the origin O' to both sides of the eUtptical sphere shadow. The distance 
between the tangent points vv:as. bisected and the resulting point connected to 
O'. The inner (^') and outer {R') endpoint of nuu'or elliptical axis is found 
by intersecting this line with . the ellipse boundary. The image M* of the 3-D 
sphere^s center point M aiid the out of plane angles a and /S with their 
respective connecting lines to the. focal spot F arc included here ottly to 
iiiustrote the spatial reliationstups. 



manual image evaluation process, the eight bit uncompressed 
tifif files were imported into an image editing software (Pho- 
toshop 7.0, Adobe Software Inc., Mountain View, Califor- 
nia). 

O' (0,0,0) was defined as one of the four ccnmil. pixels on 
the even pixel matrix. Without windowing or levelihg of the 
inriage, tangents were constructed from O' to each ellipse 
from both sides, arid the resulting tangent point connection 
line-was bisected (Fig. 3). By extending the line through this 
biisection point and O'- through the elliptical shadow, the in- 
ner (5' j aiid outer (/?') point of intersection with the ellipse 
boundary was identified and their pixel coordinates recorded. 
The coordinates of the steel-wire object endpoints as as- 
sessed along their long axis were also recorded. After con- 
version into millimcters,..thc data were fed into the algorithm 
implemented in a spread sheet software (Microsoft Excel 
2000, Microsoft Corporation, Redmond, CA). True coordi- 
nates of the sphere centers (A/, , ^2*^/3) were assessed to the 
nearest 0.5 mm with a caliper. As this error propagates into : 
our calculated depth (x) errors on a one-to-one basis, our 
experimental data include an intrinsic measurement error of 
±0.5 mm. Differences between calculated and true values 
were defined as absolute depth errors, Sx. The relative depth 
errors were computed from Sx/x. Intraobserver variability 
(reproducibility) of the manual landmark identification and 
marking process was calculated firom the standard deviation 
between three assessments pii one arbitrarily selected image 
of each, sphere radius, separated by a time interval of three 
days. 

III. RESULTS 

The mean relative error in the depth averaged over all 
images obtained with one reference radius was largest 
(-30.3%) for r=0.5 nun. For all other diameters, it ranged 
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Radius of Sphere (mm) 

Pig. 4. Experimental relative, depth eiiore presented in box plots, where 
each box defines the intecquactile distance (IQR) between the 25% and 75% 
quartile (median: bold horizxmtal. line)/' T-sluped. whiskers nuul^ values 
(ouUiers) lying «1.5 IQR beyond the quartiles; extreme values exceeding 
these Umtts are indicated by asterisks: The plots Ulusu-atc a trend of the 
underestimation of truth decreasing together with the^box heights with a 
laiger sphere radius. 



between -6.5%{r=1.5 mm) and 2.8% (r==2.5 ram). The 
95% Confidence Intervals (95% CI) were continuously de- 
creasing with increasing reference sphere radius. The results 
show a clear underestimation of depth, i.e., most of the x 
coordinates were short of the truth (Table 1; Fig. 4). Accuracy 
in X coordinates was clearly correlated with the amount of 
scatter equivalent (Spearman Rho =0.490,r<0 001; Fig. 5). 
with the absolute errors increasing with the latter. The iii|e2m 
error in object length v^^as largest (11.3%) for the smallest 
reference sphere radius (r«0.5 mm), and relatively similar 
(dkl%) for all remaining diameters. The 95% CFs of this 
error followed the same pattern as for the x coordinates 
(Table I; Fig. 6). The^ standard: deviation between single 
assessments (reproducibility) performed by oiie observer 
(F.R.) ranged between 0% and 8% of the true depth value. 

IV. DISCUSSION 

The algorithm inferes the 3-D locations of spherical ref- 
erence bodies from their elliptically distorted .,sh:^6ws. 
Spheres are frequently used in radiolbgy,^*^*^*^^'*''?-' e.g., in 
radiostercomctty (RSA);^^"^^ 3rD reconstruction from 2-D 




1.5 4.5 7.5 

Wax Plate Thickness (nim) 

FlO. 5. Box plots (for details see Fig. 4) of experimental absolute depth 
errors within the set of 0.8 s exposure time plotted against the scatter 
equivalent thickness. Except for the maximum wax plate thickness 
(7.5 mm), ogain a trend of the underestimation of troth is evident. 
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Radius of Sphere (mm) 

Fto. 6. Box plots (for details see Fig. 4) of length assessment accuracy, 
which clearly increased with the sphere radius. 



views,^'^*^ or local magnification assessment in dental 
intra-oral.'^*''*^'*^^ panoramic,* or cephalomelric^^ radiogra- 
phy. All these techniques, however, only consider 



magnification 



9-11.20.21,27^8 



or the 2-D position* *'^*^^'* of the 



marker image, yet not its elliptical distortion. This geometric 
information together with the knowledge of one distance in 
space (the sourcc-to-rcceptor distance), however, is sufficient 
to compute the sphere*s center point*s 3-D coordinates 
within the projection geometry frpm a single projection. Un- 
like the method described by Hoffmann and Esthappan,^^ the 
. aljgpntlun computes the correspondence between object and 
image points firom large image features (the sphere shadows) 
by means of geometric a priori knowledge. Hence, the iden- 
tification of image landniark points should be substantially 
easier than : in other methods.' ''•^^ Three noncollinear 
spheres attached to the pbject under smdy provide three 
known ^object" points in spacCj thereby constituting its pos- 
sible six degrees of freedom.^* If not known a priori, the 
point of incidence of the central x ray on the receptor plane 
may be reconstructed, from the point of closest approach ^be- 
tween the extended lines through at least two major elliptical 
shadow axes. Interestingly, another device especially con- 
structed for an accurate location of the point of the central 
x-ray incidence is aliso based on the divergent projection vec- 
tors on a flat image plane.'^ Once the imaging geometry is 
.known, the 3-D reconstruction of the information inherent in 
at least two 2-D views obtained uiider different geometries is 
a straightforward task.^ * Hence, the algorilhro may be used 
for image registration prior to the reconstruction process. 

If an object (distance) is coplanar with the sphere triplet, 
its true dimensions may be obtained from simple backprojec- 
tion. This may be used, if the spatial orientation of the de- 
sired dimension within the human . body can be roughly esti- 
mated a priori. In this case. Eqs. (9>-(l 1) provide the true 
3-D coordinates of any point identifiable in the image. A 
possible extension would be to attach several reference 
spheres roughly, following the object shape, allowing for the 
formation of multiple triplets coplanar to specific parts of the 
object. This would make sense for the evaluation of larger 
structures (for instance, skull or chest radiographs), where 
more than one internal structure is under examination. 
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Table I. The relative error tn the calculated reference sphere depth (x) coordinates and resulting errors in the 
calculated 6bjeet length for the four reference sphere radii. Bold values highlight the smallest errors in the 
respective row. CI: Confidence Interval. 


Sphere mdius (mm) 


0.5 


1.0 


1.5 


2.5 


Jf-cooTdinatc % error 
Mean .value 
95%-CI 

Object-length % enor 
Mean value 
95% -CI 


-30,3 
-139.4-78.8 

11.3 
-16.7-39.3 


-4.9 
-35.7-25.9 

0.6 
-8.1-9.3 


-6.5 
-34.2-21.2 

0.4 

-5.4-6.2 


2.8 
-11.6-17.2 

-1,0 
-4.5-2.5 



Operational problems include occlusion by the. sphere 
shadows, particularly in small aiea deriial- radiographs,, and 
difficulties in- attaching the spheres to small-size objects: such 
as teeth. Ho\yever, the latter has been demonsti^ed v 
elsewhere* ^'^^ Apart from geometric unsharpn^s and other 
sources of noise, geomedric eiirbrs y^ll affect tKe method -s 
accuracy in locatmg the spheres in three dimensions (for a 
detailed analysis see the internet EPAPS version of Oiis docu- 
ment). Small errors mxp and in the location of O' will have 
a negligible influence. However, an accurate identification of 
the image points R' and 5' is crucial, since they determine 
the angle a between the tangent projection lines, into which 
the algorithm fits the known sphere size. The simpliflied case 
of an orthogonal projection 08=0), where xy^Xfr-FM and 
Eq. y'\e\d^ Xf/=^Xf-rH\n{oil2)\ also clearly illustrates the 
influence of the sphere radius r. One pixel error propagates 
into a depth error rang^ig between 14%(r=2,0 mm) and 
40%(r=0.5 mm) for a typical . dental cpiifiguration {xp 
= 200 mm; pixel size: 0,039 mm) vs 3?%(r=5.0 nim) to 
75 % (r=2.0 mm) for a typical medical radiographic situation 
(jc/r=^1000 mm; pixel size: 0.130 mm). A small pixel size, 
i.e., a high resolution, reduces errors caused by the discrete 
image nature. Subpixel accuracy , in the assessment of the 
critical major elliptical axis dimension should be ainied for. 
The experiment only perfonnried to prove the principle feasi- 
bility of the method yielded depth: coordinates sh ort of tnitfa, 
i.e.. an underestiinatipn of the nisypr elU^ akis |i<?viS" | most 
probably due to low4ocal cbnttast atithe shadow's boiiiid^- 
Despite an expected size pYercstfimatipn for .s™ 
due to the weU-known point spread effect,?^ pixels were .ob- 
viously incorrectly manual ly as signed to background, al- 
though in fact: belonging xo R*S\ 

Future work will focus on a software implementation of 
the algorithm to enhance accuracy. Possible solutions include 
a generalized Hough Transform^^ for automated sphere 
shadow detection in combination with a subpixel fitting pro- 
cedure similar to that described in Ref 32. Geometric un- 
sharpness could be addressed by deconvolution with a math- 
ematically modeled inverse point spread function. Also, the 
ellipse main axis dimension could be corrected by weighting 
of the gray values* first and second derivative alotig the scan 
line R'S'?^ Once an acceptable limit of accuracy is 
achieved, our research will be directed toward 3-D recon- 



structibn from 2-D views. For this purpose, the method will 
be used for simplified registration of the 2-D images in a 
niodified Metz-Fencil technique.^ 
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APPENDIX: DERIVATION OF Ea (9) 

From the general equation of a plane in three dimensions 
[Eq. (7)] and the line intercept form equations: 

y = fcr+fl (Al) 

and 

z«fcc+*, (A2) 
we obtain the x coordinate of . the point of intersection 
TX^r.J^r.^r) between the plane and the line by substituting 
and z. from (Al) and (A2) in Eq. (7) and solving for x: 

^ A^^A2k-i'Aih 

Setting jc=0,:Eqs. (A2) and (A3) yield and z^b. Insert- 
ing .the parameters h and it in Eq. (A3) as obtained from Eqs. 
(10) and (1 1), yields the final fonn of Eq. (9): 

A-iyx*-^ A^iZr' +^ 



(A4) 
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